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In the Specification: 

Please replace each paragraph starting on the page and line as indicated with the 
following: 

Pagd/^^e 21i 

FIG. 1 is a schematic description of a cross ssection of the human retina, 
showing that the retina consists of five layers of cells, primary receptors 4, bipolar 
cells 8, amacrine cells 10, and retinal ganglion cells 12. Receptors 4 shown here are 
cone and rod cells together. The rods respond to ligght intensity, rather than color. 
The mechanism of adaptation to changing intensity modeled by Dahari and Spitzer 
operates in achromatic retinal cells 12. The receptive field (RF) [RF] of both 
chromatic and achromatic ganglion cells 12 includes both on-center and off-center 
receptors 14 of the center RF area, and receptors 16 of the surround RF area. 
Her eafter, "center" and "surround" may be written as "cen" and "smd" respectively. 

P^oelO. line 17: . 

There are six (or alternatively eight, including the yellow on- and off-center 
cells) main groups of retinal ganglion cells 12 involved in color perception, 
corresponding to the three kinds of cone cells that respond to color analogously to the 
response of receptors (rod cells) 4 to intensity. An image is first processed in three 
most common color-coded channels in the retina (L'^M", M^L' and S"^(L'^M)') to three 
activation-level maps of on-center opponent Parvo long-m edium-short wavelength (P- 
LMS) [P-LMS] cells through different retinal layers originating in (receptors) cones 4 
(FIG. I). The off-center opponent cell types L'M", M"L* and S^L^M) " have similar 
RF structure but with opposite signs. The input to the cones level is the spectral 
composition of the light reaching the retina, when illumination falls on a surface of 
the objects and is reflected from it. The field of view is mapped by the three types of 
cones, L, M and S. The quantum catches of the three cone types, Lpjgment, Mpigment and 
Spigment, are typically expressed by an integration of the inner triple product of the 
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spectral sensitivity of each pigment type with the reflectance properties and with the 
spectral and intensity absorption properties of the surface, at any specific location in 
] ! the image (Wyszeki & Styles). The responses of the three cone-types, Leone, Mcone 
^ and Scone, normalized separately to a range of 0-1, are typically expressed by a Naka- 
Rushton equation as a function of their inputs Lpigmem, Mpigmem and Spjgment (Dahari 
and Spitzer; Wyszeki & Styles). The input red, green, and blue intensity values are 
treated as though they are the responses of red, green and blue cone cells. .. . 



Pagi 



M^ine 7: 



The second stage of the forward transformation deals with the simulation of 
double-opponent responses that emulate the action of (cortical) "double-opponent" 

3 color coded Parvo double opnonent LMS fPdo-LMSs) [Pdo-LMSs] cells, and 
^ includes an additional remote adaptation. The color-coded double-opponent cells are 
adapted (corrected) by a remote adaptation in a manner similar to a mechanism based 
on psychophysical findings, as shown in Singer & D'Zmura, Vision Research, vol. 34, 
pp. 3111-3126, 1994. Adaptation is also explained (but not in relation to color and do- 
remote area) in Dahari and Spitzer. Note that "remote" adaptation refers to the effect 
of regions peripheral to the "classical" RF. 




Page* 13': line 16: 



The mathematical simulation of the physiological processes related to color 
contrast perception proceeds as follows: the spatial response profiles of the two sub- 
regions of each (on- and off-) ganglion cell RF Parvo response ganglion cells (P- 
RGC), the "center" and the "surround" regions, are preferably expressed by a 
Difference of Gaussians (DOG). As in Dahari and Spitzer, the first step of the present 
invention is the transformation of each input image to an "output function" or 
^ "spectral response function" G. For each color (red, green and blue), there is a center 

spectral response function Gc (also written as Lcen, Mcen and Seen) and a surround 
spectral response function Gs (also written as Lsrnd, Msmd and Ssmd)- Each of the three 
center response functions Gc is computed by convolving the corresponding color 
image with a center local spatial Gaussian filter fc, as shown in equations 2 and 3 of 
Dahari and Spitzer, thereby producing a center smoothed image whose pixel values 



3 



S/N: 09/777,950 





29 Apri 1-2003 
Att. Dkt. 27/181 



Art Unit: 2625 



page 4 of 1 1 



are the required GcS. Similarly, each of the three surround response functions Gs is 
computed by convolving the corresponding color image with an opposite sign 
surround local spatial Gaussian filter fs, thereby producing a surround smoothed 
image whose pixel values are the required GsS, Typical values of the radii p of the 
Gaussian filters are 0.5 pixels for the center filter and 0.5-9 pixels for the surround 
filter. 

PA^e24. line 1: , 

In order to perform contrast on real images it is necessary to inversely 
transform the adapted (corrected) double-opponent responses (at any location in the 
image) into a, perceived color. Several different inverse functions, based on different 
assumptions, can be used for this purpose in preferred embodiments of the present 
invention. The calculated perceived color contrast is the color contrast that would 
stimulate the triplet of Parvo double opponent after ad aptation (Pho-^) [Pdo-a] a LMS 
cells to the same responses, with a uniform non-contrast surface present in their 
"remote" areas. The rationale that has lead to this definition is the assumption that the 
visual system interprets color contrast of an object in a way similar to the one 
described above. The following equations describe the emulation of the three main 
double-opponent and opponent color-coded cells, but alternatively all these equations 
can be applied also to the yellow double-opponent and opponent RFs. In general, the 
steps in the inverse transformation preferably include: using the adapted double- 
opponent responses to calculate a new value of a double-opponent response ("new 
do-response"), assigning the new value to the double-opponent center ("new do-center 
response'*) thus obtaining a new value of an opponent response ("new opponent 
response"), and assigning this new opponent response value to get a new opponent 
center response .thus returning to eg. L which now has "new center" values. 



color image is connected to the color contrast adaptation described above, in the sense 
that it uses a similar physiological basis for the algorithm and a similar sequence of 
steps, except that in one alternative case ("case 1 ") the adaptation is performed on the 




Intensity adaptation or contrast adjustment on the achromatic information in a 
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center and surround areas of the opponent RFs, and not on the double-opponent RFs 
as in the color contrast case. Another difference vs the color contrast adaptation is that 
in the intensity contrast adaptation of case 1, an inverse transformation may not be 
(|^ applied. Intensity adaptation can be applied independently on achromatic images 
(non-color or black-and-white pictures) or on the intensity domain of color images. 
For intensity adaptation, the method emulates two types of ganglion cells of the 
magno pathway, the on-center and nf f-center medium wavelenp th opponent (M- 
opponenf) [ M-opponentl cells. In a preferred embodiment of the intensity adaptation, 
the remote adaptation is applied to the center sub-region of opponent RFs (or 
alternatively to the center and surround sub-regions of the opponent RFs), unlike in 
color contrast, in which remote adaptation is applied to double-opponent RFs. 
Alternatively, achromatic contrast adjustment can be performed on double-opponent 
RFs too. 



Pape 2| a/rine 13: 



where the RF region is preferably circular, fc is preferably a Gaussian decaying 
spatial-weight function, analogous to that in eq. 4. The "surround" response, which 
now represents the surround sub-region of the RF of [an] a Magno response ganglion 
cell fM-RGC) [ M-RGC], is defined as 




Paj^30. line 1: 



where fs is also a Gaussian spatial-weight function with asign opposite to that of fc. As 
in the color contrast case, the total weight of fo is 1, while the total weight of fs is 
l/rcenter/surround ratio^ H/CSR) f l/CSl. The steps above are repeated for off-center 
type cells. 



Please make the additional following changes, as suggested by the Examiner, in each 
paragraph starting on the page and line as indicated : 
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The first stage of the forward transformation is performed analogously to the 
top path of FIG, la. taken from [FIG. 1 ofl Dahari and Spitzer (notwi standing the fact 
that Dahari and Spitzer's paper does not contain color information), to provide, at each 
pixel, an opponent red center response (the pixel value of a "response function"), an 
opponent green center response, and an opponent blue center response. Optionally, a 
fourth, yellow opponent center response, which is calculated by simply adding the red 
and green responses, is obtained with the same forward transformation. In what 
follows, the fourth color (yellow) is optional in all transformations at the level of 
opponent and double-opponent cells. A similar transformation is used to obtain red, 
green and blue surround responses. These surround responses are subtracted from the 
center responses, in emulation of the action of "on-center" and "off-center" retinal 
ganglion cells, as follows: the green surround response is subtracted from the red 
center response; the red surround response is subtracted from the green center 
response; and the red surround response and the green surround response are averaged 
to form a yellow surround response, which is subtracted from the blue center 
response. Optionally, the blue surround response is subtracted from the yellow center 
response. ■ - — ' ' 



R*op(i^), R*op(m+), R*op(s+) (and optionally a fourth "yellow" R*op(Hm)+) are each both a 
spatial and a spectral filtered response, similar to Rop(i+> Rop(m+), Rop(s+) (and optionally 
the fourth "yellow" Rop((Hm>f) but related to different areas than the corresponding 
RopS. Rop and R*op and similar pairs below are also referred to herein as, respectively, 
"non-filtered" and "filtered" responses (in spite of tlie fact that the "non-filtered" 
values are obtained of course with operations involving filters). The relationship 
between non-filtered and filtered opponent responses is illustrated schematically in 
FIG. 4, where an exemplary (non-filtered) on-center opponent response Rop(xi, yO is 
calculated for pixel 62 centered at (x,, yi) using eq. 1, i.e. subtracting the value at 
(xi, yi) of a surround response function Gs(xi, yO from a center response function 
Gc(xi, yi). As explained above, both Gs(xi, yO and Gc(xi, y.i) are functions obtained 
by convolutions involving local filters. A filtered response Rop*(xi, yi) to be used in 
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^ eq. 2 for the coinputation of the do-center responses for the do-cell centered at 
(X,, y,), (and thus associated with Rop(x,, y,)) is calculated in general by applying the 
) I ) local filters in the convolutions yielding Gc and Gs at different (than x,, y,) locations 
^^JjJ g4 [62]. The "filtering "operation for obtaining filtered responses thus involves 
moving the filters over the relevant area of the integration and obtaining for each pixel 
a function value. The value of R*op is in general different from that of Rop, except for 
the special case when both are calculated at the same (x, y). 




Pals 19. line 4: 



As can be understood from the subscripts "c" (center), "s" (surround), "do-c" 
(double-opponent center), "do-s" (double-opponentsurround) and "do-r" or "do- 
remote" (double-opponent remote), the various receptive fields ("RF") and spatial 
convolutions extend over different numbers of pixels. FIGS. 5 and 6 show two 
different schematic descriptions of double-opponent RFs which include the remote 
influences. FIG. 5 shows an on-center cell 100 with a center RF 102 (pixel A) and a 
surround RF 104 (pixels B). The on-center cell is surrounded by 8 [9] off-center cells 
\ cells 106, each with a center RF 108 (pixel C) and a surround RF 110 (pixels D). Cell 

I \ 100 and cells 106 are opponent cells, and together constitute a double-opponent ("do") 
cell 112, on-center cell 100 constituting a "do-center" RF, and off-center cells 106 
constituting a "do-surround" RF. Do-cell 112 has a do-remote RF 114, which includes 
on-center opponent cell RFs 116, each with a center RF 118 (pixel E) and a surround 
RF 120 (pixels F). For simplification purposes, only nine such RF 116 are shown in 
the upper right corner of do-remote RF 114, with the understanding that they are 
repeated in the rest of the (unmarked) pixels. The different marking of pixels in cell 
100 and RFs 116 indicates the different values resulting from the application of 
different spatial weight functions at those locations. In the example of FIG. 5, the 
different receptive fields extend over different areas: for both on-center and off-center 
opponent cells, the center RF includes one pixel, while the surround RF extends over 
the 8 pixels immediately surrounding each center. For double-opponent cell 112, do- 
center RF 100 extends over 9 pixels and the do-surround RF (which includes the nine 
cells 106) extends over the adjacent 72 pixels. Do-remote RF 114 extends over the 
648 pixels surrounding do-cell 112. It should be clear that the different areas over 
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which the spatial weight functions are applied in the convolutions can also partially 
overlap, or be separated by gaps. For example, cells 116 [106] whose RFs compose 
the do-remote RF and who are used in computing the do-remote signal can partially 
overlap, or can be separated by gaps. 




Pleas^dd in p. 13 a new figure caption as follows: 



FIG. 1 is a schematic cross section of a human retina; 

FIG, la is a prior art schematic block diagram of a model of ganglio n cell response; 
FIG. 2 is a schematic description of on-center and off-center opponent receptive 
fields; 



In the Figures: 



Please add the referenced FIG. 1 of Dahari and Spitzer (mentioned on Page 11, line 17 
of the specification) as FIG. la (prior art). 

In the Claims: 

1 . (Amended) A method for correcting the color contrast of a scene, the scene 
including an intensity spectrum at each of a plurality of pixels, the method comprising 
the steps of: 

a) providing a red image, a green image, and a blue image, each image having 
a pixel value at each of the plurality of pixels; 

b) computing a center red response, a center green response and a center blue 
response based on said images; 
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